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1. INTRODUCTION 

With the rapid development of multimedia data (such as audio, video, and images) and its 
transmission across unprotected channels of communication, information security has become increasingly 
important to avoid risks and data loss [1]-[3]. Image encryption is the most widely used method of 
safeguarding the privacy of image data. Because of the unique characteristics of images, like their huge data 
capacity and strong correlation among pixels [4], [5], several conventional data encryption methods, like 
advanced encryption standard (AES) data encryption standard and (DES), are not suitable for encrypting 
digital images [6]-[8]. Therefore, various approaches for encrypting image data have been proposed, the most 
powerful and widely used of which is chaotic-based cryptography, due to its sensitivity to initial state [6], 
non-linearity, randomness, ergodicity, and so on [9]-[11], which satisfies the fundamental requirements of 
cryptography. 

In 1998, Fridrich was the first to apply the chaotic system and cryptography theory to the Image 
Encryption System, which consists of two phases: confusion and diffusion, and this architecture has become 
the most common and has been adopted by many researchers [12]. The use of deoxyribonucleic acid (DNA) 
and chaotic system-based combined cryptography to achieve high levels of security, particularly for color 
images and videos, is a broad field. The DNA molecule has several advantages, including strong parallel 
processing ability, very-low power consumption, and a high data storage capacity [13]-[15]. 
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Recently, many studies have utilized chaotic systems and DNA coding to encrypt images. Cui et al., 
[16] suggested an image encryption technique built on the hyper-chaotic-Lorenz model and dynamic DNA 
encoding. This approach uses the secure-hashing-algorithm-SHA-256 technique to produce a shared key in a 
hyper-chaotic system to create a DNA sequence matrix and a dynamic S-box. Then use SCAN mode and the 
S-box to confuse the pixel position. The pixel values are scrambled using the DNA sequence array to 
generate cipher text. Simulation results show that the technique makes keyspace more sensitive and can 
withstand statistical and differential assaults. Shima and Thanikaiselvan [17] suggested a method of image 
encryption that makes use of four different algorithms and multiple keys by splitting a single image into four 
blocks and encrypting them first with DNA as a secret key, then using RSA, then DES, and finally using the 
Chebyshev algorithm, which results in a high level of security and good image quality. 

H. Amani and M. Yaghoobi [18] propose a new 4D-hyper-chaotic dynamic system and a DNA 
sequence for color image encryption. The two key steps in this strategy are: The logical arrangement of the 
pixels is changed by employing Arnold's cat map in the first stage; permuted images are then encrypted by 
employing a combination of DNA sequences with Chen's hyper-chaotic scheme, and an adaptive technique 
makes this algorithm more complex. The results show that the method that was suggested can protect against 
brute force and statistical attacks. Huang et al., [19] suggested a strategy for encrypting chaotic color images 
using DNA-coding computations and galois field arithmetic. Firstly, three modified (1D) chaotic maps are 
presented, each with a bigger key space and improved chaotic features. Second, DNA coding and 
mathematics are used to increase the cryptosystem's permutation. Finally, the numeration applied to the 
Galois field provides the effect of pixel diffusion. It shows that the encryption method suggested has a big 
impact on how well it works, and the numerical results show that it's more secure than some of the most 
recent cryptosystems. 

Because employing only DNA to encode images is insecure, researchers are combining DNA coding 
with chaotic-systems to improve more secure image encryption techniques. However, image encryption made 
up of DNA encoding and chaotic systems has certain drawbacks, like the independency of the key generation 
method from the original image, the limitations of DNA operations, and the slow encryption speed [20]. On 
the other hand, security concerns in low-dimensional chaos have their drawbacks in the small key size and 
low security. The key size has to be large enough in order to frustrate brute force attacks [21]. As a result, 
high-dimensional chaotic systems offer higher security and randomness. In general, chaotic systems include 
only one positive Lyapunov exponent, but when the chaotic systems have at least two or more positive 
Lyapunov exponents, they are said to be hyper-chaotic. Hyper-chaotic systems generally have more complex 
dynamical behaviors than ordinary chaotic systems [22]. 

To solve the problems mentioned above, a novel 4D- hyperchaotic system and DNA encoding have 
been suggested. Firstly, the pixel positions are permuted by chaotic sequences. Secondly, according to the 
concept of DNA cryptography, a number of operations like (DNA addition, DNA XOR, DNA subtraction, 
shift right, and shift left) are applied to the DNA encoding sequence to encrypt color images and achieve 
good encryption results. 


2. THEORETICAL BACKGROUND 
2.1. The novel chaotic system 

A new 4D-hyper chaotic system is constructed using a system model represented by the following 
differential equation: 


G=-ax-bwt+cyz+ze? 
dy _ _ = z 
qo ay Tex fxz-xe (1) 


dz +h 
eg zthxy 


Wp wtixzt+jyz 

Where x,y,z, and w called the states of system, teR and a,b,c,d,e,f,g,h,i,and j are positive 
parameters of the system (1) shows a chaotic attractor in a new four-dimensional chaotic system with 
parameter values of: a=3.1, b=2.1, c=15.8, d=1.1, e=16.5, f=1.5, g=2.4, h=26.6, i=5.1, and j=12.9, and the 
initial states of : x(0)=0.2, y(0)=0.4, z(0)=1.5, and w(0)=0.8. The four Lyapunov Exponents of the nonlinear 
dynamical system (1), with parameters a=3.1, b=2.1, c=15.8, d=1.1, e=16.5, f=1.5, g=2.4, h=26.6, i=5.1, and 
j=12.9 are obtained as follows: LE;= 4.05761, LE2= 0.347562, LE3=-3.94257 and LEs=-6.61896. This 
system contains chaotic features since the highest Lyapunov exponent is positive. Because LE] and LE2 are 
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positive Lyapunov exponents, and the two remainders are negative. As a result, the system is hyper-chaotic. 
Figure 1(a) displays the system's attractors in (x-y-z), Figure 1(b) in (y-x-w), and Figure 1(c) (w-z-z). 


Figure 1. Chaotic attractors, 3-D view (a) (x-y-z), (b) (y-x-w) and (c) (x-z-w) 


2.2. DNA computing 

Adelman was the first to use DNA computing in cryptography in 1994 to solve the Hamiltonian path 
problem instead of conventional cryptography. DNA is a molecule that includes the genetic information 
necessary for any living organism's growth, development, and reproduction. DNA sequences contain of four 
nucleic acid bases: A (adenine), G (guanine), C (cytosine) and thymine (T). The complementary pairing 
concept is followed by the nucleotides Adenine and Thymine, Guanine and Cytosine. While A, C, G, and T 
indicate numerals, they can also represent the decimal digits 0, 1, 2, 3. These four decimal digits are 
equivalent to the binary numbers 00, 01, 10, and 11, which means that each nucleotide can carry two bits of 
information. This coding has a total of 4!=24 different encoding methods. However, only eight types of 
encoding schemes that refer to the complimentary pairing requirement [16], [23], [24] are shown in Table 1. 
Additionally, Table 2 includes operations for DNA addition, subtraction, and exclusive-or. 


Table 1. DNA rules [23] 


Rules 1 2 3 4 5 6 7 8 

0 A A Cc Cc G G ji E 

1 G Cc T A T A G C 

2 Cc G A T A T C G 

3 T T G G C C A A 

Table 2. The operation of DNA sequences [23] 
Addition-DNA Subtraction-DNA Exclusive-or 

(Add) A C G T (Sub) A C G T (Xor) A C G T 
A T A C G A C A T G A A C G T 
C A Cc G T C G C A T C C A T G 
G C G T A G T G Cc A G G T A Cc 
T G T A C T A T G C T T G A C 


3. METHODOLOGY 


The proposed method contains two main stages encryption stage (forward stage) and decryption 
stage (backward stage) as explained in the following paragraphs: 


3.1. Encryption stage 

The encryption stage passes several steps to produce an encrypted image, starting by splitting the 
color image into RGB bands (red, green, and blue). At the same time, the initial conditions and parameters 
are entered into the 4D-chaotic system to generate four chaotic sequences. These sequences are used in the 
permutation step for each color band. Then DNA encoding is applied to the permutation image produced 
from the previous step to use in the round function step, which includes operations such as addition, 
subtraction, and exclusive-or. The fourth chaotic sequence is used for dynamic shifting in the left and right 
directions. Three rounds are applied to find the final encrypted image as shown in Figure 2. 
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Figure 2. The general framework of encryption stage 


3.1.1. Generating 4D chaotic sequence 

Depending on the 4D chaotic system, which is generated by entering the parameters and initial 
conditions, the suggested hyper-chaotic iterates to produce four chaotic- sequences (xi, yi, zi, and wi) of real 
numbers, which are converted into three vectors as keys (key1, key2, and key3). The size of each vector is 
equivalent to the size of the original image dimensions (h*w). Table 3 shows a sample of chaotic sequences 
of real numbers that have been processed to be converted into digits of hexadecimal. 


3.1.2. Permutation step 

At this stage, the split RGB bands are reshaped into vectors (VR, VG, and VB) and are scrambled 
based on the chaotic sequences (xi, yi, zi). Each one of the three chaotic vectors will be sorted in ascending 
order and generate the index for each chaotic sequence (key index, kKey2index, Key3index) in order to permute the 
position of each value in the image bands. The first sequence (key 1 index) permutes the red vector (VR1), the 
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second sequence (key2index) permutes the green vector (VG1), and the third sequence (key3 index) permutes the 
blue vector (VB1). 


Table 3. Samples of 4D chaotic sequences 


Generated numbers Moving floating point Rounding numbers Hexadecimal form 
0.201437280655153 2014372806.5515 2014372807 7810E3C7 
0.202875019921153 2028750199.2115 2028750199 78EC4577 
0.2043 13227427431 2043 132274.2743 2043132274 79C7B972 
0.205751912809594 2057519128.0959 2057519128 7AA34018 
0.207191085709668 2071910857.0967 2071910857 TBTED9C9 
0.208630755776090 2086307557.7609 2086307558 TCSA86E6 


3.1.3. DNA encoding 

The DNA coding is applied to the result of permutation produced from the previous step by 
combining three vectors into one vector and converting each value in the vector into binary form, with every 
two bit converted into one of the DNA symbols A, C, G, and T based on Table 1. The total vector will be a 
strand of DNA sequence (Vrgbpna) in order to apply it to the DNA operations in the next steps. 


3.1.4. Round function of DNA operation 

In this step, the strand of DNA sequence (Vrgbpna) is used in DNA operations (addition, 
subtraction, and exclusive-or). Firstly, a DNA-addition operation is applied between the DNA sequence and 
the first key (key1pna) (depending on Table 2). Then, we apply the left shift to the addition result depending 
on the fourth chaotic sequence (wi). Secondly, we performed the subtraction operation between the result 
from the left shift and the second key (key2pna) (depending on Table 2), and then applied the right shift to 
the result of subtraction depending on the fourth chaotic sequence (wi). Finally, implement the Exclusive-Or 
(X-or) operation between the result right shift from the previous step and the third key vector (key3 pna). 
These steps are repeated for N rounds to obtain the final encryption image. Algorithm | and Figure 2 display 
the steps of encryption in detail. 


Encryption Algorithm (1) 
Input: original image (IM) of size hxwx3, initial conditions(x0, y0, z0, and w0) and 
parameters(a, b, c, d, e, f, g, h, I, j) 

Output: encrypted image (En) 

Begin 

step1: read (IM) image 

step2: split (IM) image into three color bands R, G, and B 

step3: assign h €&height of IM 
assign w€&width of IM 
Calculate S€ hxw 

step4: Reshape bands R, G, and B into three vectors VR, VG, and VB 

step5: Iterate suggested chaotic system (1) to generate four chaotic sequences (xi, yi, zi, 
and wi) to get(keyl, key2, key3), size of each them = S, 

step6: sort sequence xi ascending and find the position (xi) and store its position in 
index sequence(key 1 index) 
sort sequence yi ascending and find the position (yi) and store its position in 
index sequence (key? index) 
sort sequence zi ascending and find the position (zi) and store its position in 
index sequence (key3 index) 

step7: permutated VR by (Key index) to produce VR1 
Permutated VG by (key2index) to produce VG1 
Permutated VB by (key3index) to produce VB1 

step8: combine three vectors VR1, VG1, and VB1 into one vector Vrgb 

step9: Encode the vector Vrgb into DNA sequence Vrgbpna 

step10: Encode the vector key1 into DNA sequence key1 pna 
Encode the vector key2 into DNA sequence key2pna 
Encode the vector key3 into DNA sequence key3pna 

step11: Round Function: 
-apply addition operation between Vrgbpna and key1pna to produce T1 
-apply left shift T1 depending on wi sequence (wi mod sk1) to obtain TSleft 


Chaotic-DNA system for efficient image encryption ... (Huda Rashid Shakir) 


2650 O 


-apply subtraction operation between TSleft and Skey2 to produce T2 
-apply right shift T2 depending on wi sequence (wi mod sk2) to obtain TSright 
-apply Exclusive-Or operation between TSright and Skey3 to produce T3 
step12: repeat step11 for N round function 
step13: decode DNA sequence T3 to produce Tfinai 


step14: split (Trina!) into three vectors and reshape each vector into matrices (MR, MG, and 


MB) 
step16: concatenate (MR, MG, and MB) to obtained 3D matrix encrypted image (En) 
End Algorithm 


3.2. Decryption stage 


ISSN: 2302-9285 


In the decryption stage, all steps are applied to the encrypted image in reverse order. The encrypted 
image is encoded to a DNA sequence and DNA operations are applied in reverse order, which includes: (X- 
or, right shift, addition, left shift, and subtraction, respectively), and the number of rounds function is equal to 
the number of rounds in the encryption stage. The extracted image is then obtained by applying the inverse 
permutation to the result of the previous step. Algorithm 2 and Figure 3 show the decryption steps in detail. 


Decryption Algorithm (2) 
Input: Encryption image (En) of size hxwx3, initial conditions(x0, y0, z0, and w0) and 
parameters (a, b, c, d, e, f, g, h, I, j) 
Output: Extracted image (EI) 
Begin 
step1: read (En) image 
step2: split (En) image into three color bands R, G, and B 
step3: assign hCheight of En 
assign w€ width of En 
Calculate S€ hxw 
step4: Reshape three color bands into three vectors IVr, IVg, and IVg 
step5: Iterate suggested chaotic system (1) to generate three chaotic sequences (xi, yi, 
zi, and wi) to get(key1, key2, key3), size of each them > S, 
step6: combine three vectors IVg, IVc, and IVg into one vector [Vrgb 
step7: Encode the vector [Vrgb in DNA sequence IVrgbpna 
step8: Encode the vector key! in DNA sequence key 1 pna 
Encode the vector key2 in DNA sequence key2pna 
Encode the vector key3 in DNA sequence key3pna 
step9: Inverse Round function: 
-apply Exclusive-Or operation between IVrgbpna and key3pna to produce Tx1 
-apply Right shift on Tx1 depending on wi sequence (wi mod sk2) to obtain TXrignt 
-apply subtraction operation between Tx;ign and key2pna to produce Tx2 
-apply Left shift on Tx2 depending on wi sequence (wi mod sk1) to obtain Txtert 
-apply addition operation between Txef and key! pna to produce Tx3 
step10: repeat step9 for N round function 
step11: decode DNA sequence T3 to produce Terina 
step12: split (Tefinai) into three vectors IVr, IVa, and IVg 
step13: sort sequence xi ascending and find the position (xi) and store its position in 
index sequence(key lindex) 
sort sequence yi ascending and find the position (yi) and store its position in 
index sequence (key2index) 
sort sequence zi ascending and find the position (zi) and store its position in 
index sequence (key3index) 
step14: permutated the vector [Vp by key index to produce IVg: 
Permutated the vector [Vg by key2 index to produce IVa 
Permutated the vector IVg by key3 index to produce IVg: 
step15: reshape each vector (IVri, [Vci, IVB1) into matrices (MeR, MeG, and MeB) 
step16: concatenate (MeR, MeG, and MeB) to obtained 3D matrix extracted image (EI) 
End Algorithm 
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Figure 3. The general framework of decryption stage 


4. RESULTS AND DISCUSSION 

Several tests are applied to the proposed method to evaluate the efficiency of the encryption 
suggested algorithm. These tests include key space, histograms, correlation analysis. A set of color images is 
used for testing and evaluation of results. These tests are illustrated as follows: 


4.1. Key space analysis 

The key space should be large enough to thwart brute force attacks. An encryption method can use 
all possible and distinct key values as long as the minimum possible encryption key size is 10128 bits to 
withstand brute force attacks [23], [25]. In the suggested method, the key space is equal to 2°’, which means 
it is resistant to brute-force attack. Table 4 compares the proposed method's key space results with those of 
related works. 


Chaotic-DNA system for efficient image encryption ... (Huda Rashid Shakir) 


2652 O ISSN: 2302-9285 


Table 4. Comparative of key space 


Method Key space value 
Proposed method ge 
Ref.[16] 3.4028 * 10” 
Ref.[18] Dua 
Ref.[19] 222 


4.2. Histogram analysis 

The histogram shows the visual distribution of pixel values in an image, which is an 
important aspect to reflect if the algorithm is resistant to statistical analysis. The histogram of the 
cipher image must be uniform and varied from the corresponding histogram of the original image, 
which consequently provides no clue for statistical analysis of the encrypted image [26], [27]. 
Figure 4 shows the histogram tested images before and after encryption. 


Figure 4. The histogram analysis of proposed method 


4.3. Correlation analysis 

In general, an image is made up of many pixels. Because adjacent pixels are usually quite similar, 
they have a high correlation with one another. An ideal image encryption technique would disrupt the 
intrinsic association between neighboring pixels, making it impossible for an opponent to guess the original 
image using this information. Normally, the following equation is used to determine the correlation between 
vertical, horizontal, and diagonal pixels [28]: 


N Yjea(xpxyj)—Vjea xX 


CC = 
[wets F- CFL DN EI 9F- OL 7) 


(2) 


The pixel intensity values of two neighboring pixels are x and y, and the number of pixels is denoted 
by N. Table 5 shows the correlation test for the plain and encrypted image. Figure 5 shows the correlation 
test of the suggested method in the (horizontal, vertical, and diagonal) directions. Also, Table 6 displays the 
comparison results of the correlation test of Lena image with other methods. 


Table 5. Correlation test for plain and encrypted image 


Horizontal Correlation Vertical correlation Diagonal correlation 
Plain image Encrypted image Plain image Encrypted image Plain image Encrypted image 
0.89285 0.00642 0.89285 0.00642 0.84829 0.00110 
0.92566 0.00320 0.92566 0.00320 0.86998 -0.00024 
0.98816 0.00423 0.98816 0.00423 0.97477 -0.00504 
0.96438 0.00305 0.96438 0.00305 0.94832 -0.00191 
0.92188 -0.00211 0.92188 -0.00211 0.85933 -0.00360 
0.96362 -0.00126 0.96362 -0.00126 0.92952 -0.00023 
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Figure 5. The correlation test of proposed method (horizontal, vertical and diagonal) 


Table 6. Comparison results of correlation test lena image with other methods 


Lena Encrypted image Horizontal Vertical Diagonal 
Proposed method 0.00305 0.00305 -0.00191 
[16] -0.0046 0.0045 -0.0074 
[18] 0.00603 0.00476 0.0041 
[19] - 0.0034 0.0021 - 0.0003 


4.4. Information Entropy Analysis 

Shannon proposed information entropy in 1948, using the thermodynamic idea of entropy to identify 
the relation among probability and information redundancy in mathematical terms. For a perfect random 
image, the information entropy value is 8 [29]-[31]. In (3) is the formula for calculating the information 
entropy of a random information sequence X: 


Entropy E(X) = — Yizi p(X)log p(X) (3) 
Table 7 shows the values of the entropy test of images. The results indicate that the information 


entropy of the encryption images is near to 8, which means that the system is capable of resisting entropy 
attacks. Table 8 displays the comparison results of the entropy test of the Lena image with other methods. 


Table 7. The entropy test of images 


Image Palace Baboon Vegetables Lena Monuments Birds 
Plain image 7.68523 7.69805 7.43719 7.75155 7.54016 6.97018 
Encryption image 7.99857 7.99656 7.99702 7.99987 7.99727 7.99740 


Table 8. Comparison results of entropy test Lena image with other methods 


Lena encrypted image Entropy value 
Proposed method 7.99987 
Ref.[16] 7.9896 
Ref.[18] 7.9983 
Ref.[19] 7.9971 
Ref.[31] 7.929 


4.5. Differential attack analysis 

The unified average changing intensity (UACI) and the number of pixels change rate (NPCR) are 
two differential attack metrics used to determine the susceptibility of the original data to minor modifications. 
Assume the cipher images are C and C' before and after modifying a single pixel in a plain image [32]-[34]. 
The following is the formula: 
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NPCR = aren x 100% (4) 
ee: IQ G@)-C2G/)| 
VACI = za [z 255 | ~ 100% ©) 


W denotes the image's width, while H denotes the image's height; C1 and C2 represent the encrypted 
images before and after changing one pixel in the original image. Table 9 shows the results of the NPCR and 
UCAI tests. In addition, Table 10 displays comparison results of the NPCR as well as UACI tests of Lena 
image with other methods. 


Table 9. The NPCR and UCAI test 


Image NPCR UCAI 
Palace 99.61090088 31.24479406 
Baboon 99.61242676 29.28272123 
vegetables 99.61547852 33.6586567 
Lena 99.60699544 33.36560258 
monuments 99.62005615 31.6231119 
birds 99.60327148 35.36180085 


Table 10. Comparison results of The NPCR and UACTI test of Lena image with other methods 


Lena encrypted image NPCR UCAI 
Proposed method 99.60699544 33.36560258 
Ref.[16] 99.60 28.71 
Ref.[17] 98.6 33.1 
Ref.[18] 99.53 34.26 
Ref.[19] 99.6108 33.4642 
Ref.[31] 99.65 32.4966 


4.6. MSE and PSNR ratio analysis 

As a common criterion for image encryption techniques, the cipher-image should be significantly 
different from the original image. The difference between plain and cipher images can be assessed using two 
criteria: the mean square error (MSE) and the peak signal-to-noise ratio (PSNR). Table 11 shows the MSE 
and PSNR tests and they can be calculated using (6) and (7) [24]. 


PSNR = 20 logio(#38,)aB (6) 


MSE = wn Žij (po (i-j) — pr (i.j))? (7) 


Table 11. The MSE and PSNR test 


Image Palace Baboon vegetables Lena Monuments Birds 
MSE 9505.137 8236.518 11357.88 8914.025 8410.986 15437.31 
PSNR 3.607009 4.242124 2.89338 3.891629 4.073957 1.440728 


5. CONCLUSION 

A novel image encryption approach is introduced using a 4D-hyper chaotic system and DNA 
encoding to achieve a high level of image encryption security. Firstly, chaotic sequences generated by a new 
4D-hyperchaotic system are used to permute the pixel positions. Secondly, the pixel values are changed 
according to DNA operations (like DNA addition, DNA XOR, DNA subtraction, shift right, and shift left) to 
achieve high diffusion. The algorithm's performance was evaluated using analysis measures such as 
histogram analysis, key space analysis, correlation coefficient analysis, information entropy analysis, number 
of pixels change rate (NPCR), unified average changing intensity (UACI), PSNR, and MSE. The results of 
the tests show that the algorithm is very secure, efficient, and resistant to a lot of different types of attacks. 
For encrypted images, the histogram is fairly uniform, the correlation values between adjacent pixels are very 
small and close to zero, and the entropy is close to the ideal value of eight. Also, the NPCR and UACTI values 
are close to the ideal values of (99%) and (33%), respectively. Additionally, the suggested system is 
extremely sensitive to key changes and has a very wide key space equivalent to 2627 keys, which makes it 
withstand brute-force, differential, as well as statistical attacks. 
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